A nonstandard finite difference scheme for a time-fractional model of Zika virus transmission.

In this work, we investigate the transmission dynamics of the Zika virus, considering both a compartmental model involving humans and mosquitoes and an extended model that introduces a non-human primate (monkey) as a second reservoir host. The novelty of our approach lies in the later generalization of the model using a fractional time derivative. The significance of this study is underscored by its contribution to understanding the complex dynamics of Zika virus transmission. Unlike previous studies, we incorporate a non-human primate reservoir host into the model, providing a more comprehensive representation of the disease spread. Our results reveal the importance of utilizing a nonstandard finite difference (NSFD) scheme to simulate the disease's dynamics accurately. This NSFD scheme ensures the positivity of the solution and captures the correct asymptotic behavior, addressing a crucial limitation of standard solvers like the Runge-Kutta Fehlberg method (ode45). The numerical simulations vividly demonstrate the advantages of our approach, particularly in terms of positivity preservation, offering a more reliable depiction of Zika virus transmission dynamics. From these findings, we draw the conclusion that considering a non-human primate reservoir host and employing an NSFD scheme significantly enhances the accuracy and reliability of modeling Zika virus transmission. Researchers and policymakers can use these insights to develop more effective strategies for disease control and prevention.

A physics-informed neural network to model COVID-19 infection and hospitalization scenarios.

In this paper, we replace the standard numerical approach of estimating parameters in a mathematical model using numerical solvers for differential equations with a physics-informed neural network (PINN). This neural network requires a sequence of time instances as direct input of the network and the numbers of susceptibles, vaccinated, infected, hospitalized, and recovered individuals per time instance to learn certain parameters of the underlying model, which are used for the loss calculations. The established model is an extended susceptible-infected-recovered (SIR) model in which the transitions between disease-related population groups, called compartments, and the physical laws of epidemic transmission dynamics are expressed by a system of ordinary differential equations (ODEs). The system of ODEs and its time derivative are included in the residual loss function of the PINN in addition to the data error between the current network output and the time series data of the compartment sizes. Further, we illustrate how this PINN approach can also be used for differential equation-based models such as the proposed extended SIR model, called SVIHR model. In a validation process, we compare the performance of the PINN with results obtained with the numerical technique of non-standard finite differences (NSFD) in generating future COVID-19 scenarios based on the parameters identified by the PINN. The used training data set covers the time between the outbreak of the pandemic in Germany and the last week of the year 2021. We obtain a two-step or hybrid approach, as the PINN is then used to generate a future COVID-19 outbreak scenario describing a possibly next pandemic wave. The week at which the prediction starts is chosen in mid-April 2022.

A nonstandard finite difference scheme for the SVICDR model to predict COVID-19 dynamics.

In the context of 2019 coronavirus disease (COVID-19), considerable attention has been paid to mathematical models for predicting country- or region-specific future pandemic developments. In this work, we developed an SVICDR model that includes a susceptible, an all-or-nothing vaccinated, an infected, an intensive care, a deceased, and a recovered compartment. It is based on the susceptible-infectious-recovered (SIR) model of Kermack and McKendrick, which is based on ordinary differential equations (ODEs). The main objective is to show the impact of parameter boundary modifications on the predicted incidence rate, taking into account recent data on Germany in the pandemic, an exponential increasing vaccination rate in the considered time window and trigonometric contact and quarantine rate functions. For the numerical solution of the ODE systems a model-specific non-standard finite difference (NSFD) scheme is designed, that preserves the positivity of solutions and yields the correct asymptotic behaviour.

Enhancing the spatial resolution of hyperpolarized carbon-13 MRI of human brain metabolism using structure guidance.

PURPOSE: Dynamic nuclear polarization is an emerging imaging method that allows noninvasive investigation of tissue metabolism. However, the relatively low metabolic spatial resolution that can be achieved limits some applications, and improving this resolution could have important implications for the technique. METHODS: We propose to enhance the 3D resolution of carbon-13 magnetic resonance imaging (13 C-MRI) using the structural information provided by hydrogen-1 MRI (1 H-MRI). The proposed approach relies on variational regularization in 3D with a directional total variation regularizer, resulting in a convex optimization problem which is robust with respect to the parameters and can efficiently be solved by many standard optimization algorithms. Validation was carried out using an in silico phantom, an in vitro phantom and in vivo data from four human volunteers. RESULTS: The clinical data used in this study were upsampled by a factor of 4 in-plane and by a factor of 15 out-of-plane, thereby revealing occult information. A key finding is that 3D super-resolution shows superior performance compared to several 2D super-resolution approaches: for example, for the in silico data, the mean-squared-error was reduced by around 40% and for all data produced increased anatomical definition of the metabolic imaging. CONCLUSION: The proposed approach generates images with enhanced anatomical resolution while largely preserving the quantitative measurements of metabolism. Although the work requires clinical validation against tissue measures of metabolism, it offers great potential in the field of 13 C-MRI and could significantly improve image quality in the future.

Equivariant neural networks for inverse problems.

In recent years the use of convolutional layers to encode an inductive bias (translational equivariance) in neural networks has proven to be a very fruitful idea. The successes of this approach have motivated a line of research into incorporating other symmetries into deep learning methods, in the form of group equivariant convolutional neural networks. Much of this work has been focused on roto-translational symmetry of R d , but other examples are the scaling symmetry of R d and rotational symmetry of the sphere. In this work, we demonstrate that group equivariant convolutional operations can naturally be incorporated into learned reconstruction methods for inverse problems that are motivated by the variational regularisation approach. Indeed, if the regularisation functional is invariant under a group symmetry, the corresponding proximal operator will satisfy an equivariance property with respect to the same group symmetry. As a result of this observation, we design learned iterative methods in which the proximal operators are modelled as group equivariant convolutional neural networks. We use roto-translationally equivariant operations in the proposed methodology and apply it to the problems of low-dose computerised tomography reconstruction and subsampled magnetic resonance imaging reconstruction. The proposed methodology is demonstrated to improve the reconstruction quality of a learned reconstruction method with a little extra computational cost at training time but without any extra cost at test time.

Synergistic multi-spectral CT reconstruction with directional total variation.

This work considers synergistic multi-spectral CT reconstruction where information from all available energy channels is combined to improve the reconstruction of each individual channel. We propose to fuse these available data (represented by a single sinogram) to obtain a polyenergetic image which keeps structural information shared by the energy channels with increased signal-to-noise ratio. This new image is used as prior information during a channel-by-channel minimization process through the directional total variation. We analyse the use of directional total variation within variational regularization and iterative regularization. Our numerical results on simulated and experimental data show improvements in terms of image quality and in computational speed. This article is part of the theme issue 'Synergistic tomographic image reconstruction: part 2'.

Motion estimation and correction for simultaneous PET/MR using SIRF and CIL.

SIRF is a powerful PET/MR image reconstruction research tool for processing data and developing new algorithms. In this research, new developments to SIRF are presented, with focus on motion estimation and correction. SIRF's recent inclusion of the adjoint of the resampling operator allows gradient propagation through resampling, enabling the MCIR technique. Another enhancement enabled registering and resampling of complex images, suitable for MRI. Furthermore, SIRF's integration with the optimization library CIL enables the use of novel algorithms. Finally, SPM is now supported, in addition to NiftyReg, for registration. Results of MR and PET MCIR reconstructions are presented, using FISTA and PDHG, respectively. These demonstrate the advantages of incorporating motion correction and variational and structural priors. This article is part of the theme issue 'Synergistic tomographic image reconstruction: part 2'.

(An overview of) Synergistic reconstruction for multimodality/multichannel imaging methods.

Imaging is omnipresent in modern society with imaging devices based on a zoo of physical principles, probing a specimen across different wavelengths, energies and time. Recent years have seen a change in the imaging landscape with more and more imaging devices combining that which previously was used separately. Motivated by these hardware developments, an ever increasing set of mathematical ideas is appearing regarding how data from different imaging modalities or channels can be synergistically combined in the image reconstruction process, exploiting structural and/or functional correlations between the multiple images. Here we review these developments, give pointers to important challenges and provide an outlook as to how the field may develop in the forthcoming years. This article is part of the theme issue 'Synergistic tomographic image reconstruction: part 1'.

Learning the Sampling Pattern for MRI.

The discovery of the theory of compressed sensing brought the realisation that many inverse problems can be solved even when measurements are "incomplete". This is particularly interesting in magnetic resonance imaging (MRI), where long acquisition times can limit its use. In this work, we consider the problem of learning a sparse sampling pattern that can be used to optimally balance acquisition time versus quality of the reconstructed image. We use a supervised learning approach, making the assumption that our training data is representative enough of new data acquisitions. We demonstrate that this is indeed the case, even if the training data consists of just 7 training pairs of measurements and ground-truth images; with a training set of brain images of size 192 by 192, for instance, one of the learned patterns samples only 35% of k-space, however results in reconstructions with mean SSIM 0.914 on a test set of similar images. The proposed framework is general enough to learn arbitrary sampling patterns, including common patterns such as Cartesian, spiral and radial sampling.

Compartment models for vaccine effectiveness and non-specific effects for Tuberculosis.

In this paper, we attempt to set a framework of conditions for model-specific predictions of newly arising TB epidemics by e.g. immigration of infected persons from high prevalence countries. In addition, we address the aspect of trained immunity in our model. Using a mathematical approach of a system of ordinary differential equations which can be developed over several time-points we obtained varying infection or attack rates that led to different effects of the vaccination, depending on the setting of certain parameters and starting values in the compartments of a SEIR-model. We finally obtained different graphs of disease progression and were able to outline which upgrades and expansions our system requires in order to be exact and well adapted for predicting the course of future TB outbreaks. The model might also be beneficial in predicting non-specific effects of vaccines.

Faster PET reconstruction with non-smooth priors by randomization and preconditioning.

Uncompressed clinical data from modern positron emission tomography (PET) scanners are very large, exceeding 350 million data points (projection bins). The last decades have seen tremendous advancements in mathematical imaging tools many of which lead to non-smooth (i.e. non-differentiable) optimization problems which are much harder to solve than smooth optimization problems. Most of these tools have not been translated to clinical PET data, as the state-of-the-art algorithms for non-smooth problems do not scale well to large data. In this work, inspired by big data machine learning applications, we use advanced randomized optimization algorithms to solve the PET reconstruction problem for a very large class of non-smooth priors which includes for example total variation, total generalized variation, directional total variation and various different physical constraints. The proposed algorithm randomly uses subsets of the data and only updates the variables associated with these. While this idea often leads to divergent algorithms, we show that the proposed algorithm does indeed converge for any proper subset selection. Numerically, we show on real PET data (FDG and florbetapir) from a Siemens Biograph mMR that about ten projections and backprojections are sufficient to solve the MAP optimisation problem related to many popular non-smooth priors; thus showing that the proposed algorithm is fast enough to bring these models into routine clinical practice.

Fast Quasi-Newton Algorithms for Penalized Reconstruction in Emission Tomography and Further Improvements via Preconditioning.

This paper reports on the feasibility of using a quasi-Newton optimization algorithm, limited-memory Broyden-Fletcher-Goldfarb-Shanno with boundary constraints (L-BFGS-B), for penalized image reconstruction problems in emission tomography (ET). For further acceleration, an additional preconditioning technique based on a diagonal approximation of the Hessian was introduced. The convergence rate of L-BFGS-B and the proposed preconditioned algorithm (L-BFGS-B-PC) was evaluated with simulated data with various factors, such as the noise level, penalty type, penalty strength and background level. Data of three 18F-FDG patient acquisitions were also reconstructed. Results showed that the proposed L-BFGS-B-PC outperforms L-BFGS-B in convergence rate for all simulated conditions and the patient data. Based on these results, L-BFGS-B-PC shows promise for clinical application.

NiftyPET: a High-throughput Software Platform for High Quantitative Accuracy and Precision PET Imaging and Analysis.

We present a standalone, scalable and high-throughput software platform for PET image reconstruction and analysis. We focus on high fidelity modelling of the acquisition processes to provide high accuracy and precision quantitative imaging, especially for large axial field of view scanners. All the core routines are implemented using parallel computing available from within the Python package NiftyPET, enabling easy access, manipulation and visualisation of data at any processing stage. The pipeline of the platform starts from MR and raw PET input data and is divided into the following processing stages: (1) list-mode data processing; (2) accurate attenuation coefficient map generation; (3) detector normalisation; (4) exact forward and back projection between sinogram and image space; (5) estimation of reduced-variance random events; (6) high accuracy fully 3D estimation of scatter events; (7) voxel-based partial volume correction; (8) region- and voxel-level image analysis. We demonstrate the advantages of this platform using an amyloid brain scan where all the processing is executed from a single and uniform computational environment in Python. The high accuracy acquisition modelling is achieved through span-1 (no axial compression) ray tracing for true, random and scatter events. Furthermore, the platform offers uncertainty estimation of any image derived statistic to facilitate robust tracking of subtle physiological changes in longitudinal studies. The platform also supports the development of new reconstruction and analysis algorithms through restricting the axial field of view to any set of rings covering a region of interest and thus performing fully 3D reconstruction and corrections using real data significantly faster. All the software is available as open source with the accompanying wiki-page and test data.

PET Reconstruction With an Anatomical MRI Prior Using Parallel Level Sets.

The combination of positron emission tomography (PET) and magnetic resonance imaging (MRI) offers unique possibilities. In this paper we aim to exploit the high spatial resolution of MRI to enhance the reconstruction of simultaneously acquired PET data. We propose a new prior to incorporate structural side information into a maximum a posteriori reconstruction. The new prior combines the strengths of previously proposed priors for the same problem: it is very efficient in guiding the reconstruction at edges available from the side information and it reduces locally to edge-preserving total variation in the degenerate case when no structural information is available. In addition, this prior is segmentation-free, convex and no a priori assumptions are made on the correlation of edge directions of the PET and MRI images. We present results for a simulated brain phantom and for real data acquired by the Siemens Biograph mMR for a hardware phantom and a clinical scan. The results from simulations show that the new prior has a better trade-off between enhancing common anatomical boundaries and preserving unique features than several other priors. Moreover, it has a better mean absolute bias-to-mean standard deviation trade-off and yields reconstructions with superior relative l2-error and structural similarity index. These findings are underpinned by the real data results from a hardware phantom and a clinical patient confirming that the new prior is capable of promoting well-defined anatomical boundaries.

Vector-valued image processing by parallel level sets.

Vector-valued images such as RGB color images or multimodal medical images show a strong interchannel correlation, which is not exploited by most image processing tools. We propose a new notion of treating vector-valued images which is based on the angle between the spatial gradients of their channels. Through minimizing a cost functional that penalizes large angles, images with parallel level sets can be obtained. After formally introducing this idea and the corresponding cost functionals, we discuss their Gâteaux derivatives that lead to a diffusion-like gradient descent scheme. We illustrate the properties of this cost functional by several examples in denoising and demosaicking of RGB color images. They show that parallel level sets are a suitable concept for color image enhancement. Demosaicking with parallel level sets gives visually perfect results for low noise levels. Furthermore, the proposed functional yields sharper images than the other approaches in comparison.